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V. On a Class of Differential Equations, including those which occur in Dynamical. 
Problems.— Part I. By W. F. Donkin, M.A., F.R.S., F.R.J.S., Savilian Pro- 
fessor of Astronomy in the University of Oxford. 

Received February 23, — Read February 23, 1854. 

1 HE Analytical Theory of Dynamics, as it exists at present, is due mainly to the 
labours of Lagrange, Poisson, Sir W. R. Hamilton, and Jacobi ; whose researches 
on this subject present a series of discoveries hardly paralleled, for their elegance 
and importance, in any other branch of mathematics. 

The following investigations in the same department do not pretend to make any 
important step in advance ; though I should not of course have presumed to lay 
them before the Society, if I had not hoped they might be found to possess some 
degree of novelty and interest*. 

Of previous publications with which I am acquainted, those most nearly on the same 
subject are, Sir W. R. Hamilton's two memoirs " On a General Method in Dynamics" 
in the Philosophical Transactions; Jacobi's Memoir in the 17th vol. of Crelle's 
Journal, " Ueber die Reduction der partiellen Differential-gleichungen," &c. ; and 
M. Bertrand's " Mtimoire sur l'integration des equations difF6rentielles de la M6ca- 
nique," in Liouville's Journal (1852). The relation in which the present essay stands 
to the papers just named will be apparent to those who are acquainted with them, 
and it would be useless to attempt to make it intelligible to others. 



Oxford, Feb. 21, 1854. 



Section I. 



1. Let j?„ %, .... x n be n variables, connected by n relations with n other variables 
y!,y 2 , ...y n ; so that each variable of either set may be considered as a function of 
the variables of the other set. Suppose then 

yi == Pi(Xi, x 2 , .... x n ), 

[* It may be useful to specify the parts to which I should principally refer as containing what is, relatively 
to my own reading on the subject, new ; and in the present day it can hardly be required of any one to profess 
more than this kind of originality. These are — the theorem (3.), art. 1. The results of arts. 2 to 4. The 
formulae (19.), art. 7. The general form of the theorem (26.), art. 10. The processes and results of arts. 12 
to 14. The generalization of Sir W. Hamilton's transformation of the dynamical equations, arts. 17, 18. 
The demonstration of Poisson's theorem, arts. 21, 22. The contents of art. 25. The method of obtaining 
elliptic elements, arts. 27 to 30. The contents of arts. 34 to 36. The solution of the problem of rotation, 
Section III.] 
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this equation would become identical if x x , x m ... x n , in its second member, were 
expressed in terms of y u y%, ... y n ; hence, differentiating each side, on this hypothesis, 
first with respect to y h and then with respect to y } , we obtain 

dx x dyi'd% 3 dyi''"'da! n dyt •••••■••• v •/ 
dx 1 dyj'dx^dyj''"'dx n dy/ •-••••••{■} 

where j is any index different from i. These theorems are given by Jacobi in his 
memoir "De Determinantibus functionalibus." They are however only particular 
cases of more general theorems, which may be investigated as follows. 
If we represent by 

h J> "5 •■•• 

p, q, r, .... 

any two determinate sets of m indices each, selected out of the series 1, 2, 3, ... n, 
then the determinant formed with the m 2 differential coefficients 

dyi dy { dyj d Vj . & 

ViObjy tttt/q ihtki.n GitX/q 

possesses properties remarkably analogous to those of a simple differential coefficient. 
This analogy was pointed out by Jacobi, and has been further developed by M. Ber- 
trand in his " M^moire sur le Determinant d'un systeme de Fonctions" (Liouville's 
Journal, 1851). 

It appears to me that such functional determinants might be appropriately and 
conveniently denoted by a symbol analogous to that of a common differential coeffi- 
cient ; thus 

fe»»y*--->, ( D .) 

™\pPpy <£q} $ry • • •/ 

and I shall adopt this notation in the present paper. For example, 

d(u, v) 
d(x, y) 

would represent the determinant 

du dv du dv 
dx dy dy dx 

[The expression (D.) is not a mere arbitrary symbol, but, like a simple differential 
coefficient, is a real fraction. For if we denote by 

U\Xp) Xg, X,., ...J 

the determinant formed with the m 2 quantities 

"l**«5 ™V^qj f^l*X/yy • • • • 

t*2«*pj &<%ju g , U>2^ry • ♦ * 



^m^pi U'mpVqU'my &ry •••* 
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and attribute a corresponding meaning to 

where d x , d 2 , ... d n are symbols denoting n distinct and independent sets of variations, 
so that 

dv-^dx+^dx4- 4-^djv 

then it follows from well-known properties of determinants (as M. Bertrand has 
shown) that the complete functional determinant formed with the n 2 differential 
coefficients 

%i, dy K dy, dy„ - 

is equal to the quotient of the two determinants which I propose to denote by 

'•(i/l) 3^25 3^3) ••• Jfnji u\X l} x 2) x 3) ... x n ) } 

and moreover that the partial functional determinant formed with the m 2 terms 

%,%£ %, fa, ... &c. 

(AjJUp U/iJLfi VbSC-n iiiiAi q 

is equal to the quotient of the two partial determinants 

the differentials of y t , &c. being taken on the hypothesis that all the differentials of 
the ^-variables are =0, except those of the set x p , x q , x r , .... Thus the expression 
(D.) is a real fraction, provided its numerator and denominator be interpreted in a 
manner exactly analogous to that in which the numerator and denominator of an 
ordinary total or partial differential coefficient are interpreted.] 

This being premised, let «„ u 2 , .... u m be m functions of any or all of the functions 
y l} y 2 , ... y n (m being supposed not greater than n), so that u x , u 2 , &c. are functions 
of x„ x 2 , &c. through y u y 2 , &c. 

Let any selected sets of m indices out of the series 1,2, ... n, be denoted, for greater 
clearness, by u l} u 3 , ... « m ; (3 15 13 2 , ... (3 ro ,&c. Then the general theorem analogous to 

may be expressed as follows : — 

(the summation on the second side referring only to the indices |3, and extending to 
every combination of m out of the n numbers 1, 2, ... n). 
In like manner, the theorem analogous to 

dui dui dy x dui dy% | 

dxj dy x dxj ' dy^ dxj < "" 

<*(««,> «*,>••• ««J _ v /«*(««,> »«,,... u a j <%,, y fc , . . . ypj 

IS -J7 r — -£, 



'13 



d(x yi , x w . . . x Ym ) L% Pl , y fe , . . . ypj rf(a? Yl , a? w . . . x y J 

MDCCCLIV. L 
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These two theorems (expressed in a different notation) may be found in the memoirs 
above cited. But the following, which we shall have occasion to employ hereafter, 
has not, so far as I am aware, been explicitly stated. 

d/it' ti/ij ' 

Inasmuch as :£r.= lj ~7hJ = Q> i fc follows that the determinant represented by 

is =1 if |8 13 (3 g , ... |8 m be the same combination of indices as «,, a 2 , ... a m , but is =0 in 
every other case. (For in the first case the determinant is formed with 1, 0, 0, ... ; 
0, 1, 0, — ; 0, 0, 1, ... ; &c.j but if there be one index p^ which is not contained in 
the series «„. « 2 , &c, then one row of terms in the determinant will consist wholly of 
zeros.) 

Now considering - y 1} y 2 , &c. as functions of x u x 2 , &c, and again considering these 
latter as functions of y u y 2 , &c. given by the inverse equations, we have, by the prece- 
ding theorem, for the value of the determinant (E.) above written, the expression y„= 

\"'{yct. l , yx0 ••• yam) d(sc yi , a? y2 , ... x 7m ) J 



y \d(x yi , x y2 , ... x ym ) % ft , y A , ... ypjj 

(where «„ oj 2! ... a„; j3 15 j8 s , ... /3 m are two determinate sets of m out of the n indices, 
and the summation with respect to the indices y extends to every combination of m 
out of the ra). Consequently, 

v-= l or v»=o, (3.) 

according as the series of indices 

Pi, P25 ••• Pm 

is, or is not, the same combination as 

a i> a 25 • • • • a m- 

(I suppose, for convenience, that when the two combinations are the same, the arrange- 
ment is the same in each ; otherwise the value of y m may be — 1.) 

This is the theorem in question. If we put m=\, we obtain the equations (1.) and 
(2.) given at the beginning of this article. If we put m=n, the expression y B reduces 
itself to the product of the two determinants formed respectively with the complete 

fill: (JT • 

sets of differential coefficients -jg &c, j 2 &c, the value of which product is = 1, as is 

well known. 

As an illustration , it may be useful to exhibit the theorem in the case of m=2, 
as expressed by the common notation. Namely, 

V2 ~ *\\dg dx—dxjdxi) \dy x dy~dy ? dy x )]— 1 > 0V ~ U > ' ' • • t 4 

according as a, (S are, or are not, the same as p, q. Here a, (3 ; p, q are two deter- 
minate pairs of indices, and the summation refers to i, j, extending to every binary 
combination. 
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2. Theorem. — Retaining the suppositions made at the beginning of the last article, 
let X be a given function of x x , x 2 , .... x n ; and let us further suppose that the equations 
by which y 1} y 2 , -•■yn are determined as functions of j?„ &c, are 

_dX _dX _dX ) / 5 ^ 

(130 \ CttSOn UiS/yi 

so that &=&; 

(130 j (tOC i 

and if we transform the equations (1.), (2.), art. 1, by this condition, we obtain the 
n equations 

dx l dyi* dx 1 dy i ' r '"*'dx 1 dy { 

d]h dxj difa dx% dy», dx* 

dx s dyi ' dx^ dyC* " dx% dy { 

djh do^.dy^ dx% dy^ dxn__ 

dxi dy t 'dxi dy^ '"' dx t dy { 

If these equations be added, after multiplying them respectively by 

dyj dyj dyj 

the sum of the first members reduces itself by virtue of the equations (1.), (2.), to j- s , 

whilst the second side consists of the single term -—■• We have then 



dy~dyi 

or, in other words, j/"x„ x 2 , ... x n be found from the system of equations (5.) in terms of 
yi> y2) ••• y™ the resulting expressions are the partial differential coefficients of a certain 
function of 'y \, y 2 , ..., y n , so that the system inverse to (5.) is of the form 

» _^Y „ _dY _dY , R v 

dy x dy z ' % v ' 

The relation between X and Y is easily found as follows. The equations (5.) and (6.) 

give 

dX=y l dx l +y 2 dx 2 + ... +y n dx n 

dY=x l dy l +x 2 dy 2 + . . . +x n dy n ; 

whence, by addition, d(X-\-Y)=d(x 1 y 1 -\'X^ 2 -\- ...-\-x n y n ), 

and therefore X+Y=x 1 y l +x 2 y 2 +...+x n y„ (7.) 

(omitting the arbitrary constant, which might of course be added). 
The actual value of Y will then be 

Y=~(X) + (x l )y 1 +(x 2 )y 2 +.., + (x n )y a , ..... (8.) 
l2 
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in which the brackets indicate that x u x 2 , ... x n are to be expressed in terms of 
Vm y<a ••• Vw so that Y may be a function of the latter variables only. It is easy to 
show h posteriori that the expression (8.) verifies the equations (6.), but I pass on to 
some further considerations. (See note at the end of Section II.) 

3. Suppose the function X involves explicitly, besides the variables x u x 2 , &c, any 
other quantity p, so that the expressions (a?,), (x 2 ), &c. (or the values of x 1} x 2 , ... in 
terms of y x , y 2 , &c.) will also involve p explicitly, and we shall have 

d(X) dX dX d(xj dX % g ) 
dp dp dx x dp ' dx<i dp > "' 

_dX d(x x ) dfa) 

— dp^y^ dp ~T~y* dp ^— 

Now, differentiating the equation (8.) with respect to p (so far as it contains p 
explicitly), we obtain 

d J = _ d ^ 4 . u d M> v d M + 
dp ~ dp "■ "' dp ' " 2 dp '"' 

which the equation above written reduces simply to 

dX ,dY n . - /„ s 
!^+!p= ( 9 -) 

In the particular case in which X is a homogeneous function of x 1} x 2 , ... x n , and of 
m dimensions with respect to those variables, the equations (8.) and (9.) become 

Y=(wi-l)(X) -, 

f+(-Df =0} (10 -) 

and it is easily seen that Y is also homogeneous and of ^^ dimensions in y x , y 2 ,... y n . 

4. The theorems (8.) and (9.) are cases of more general ones which are easily proved 
in a perfectly similar way, and which I shall therefore only enunciate. If, by means 
of the equations (5.), art. 2, we express a set of n out of the 2n variables, consisting 
of r x's and n—ry's, of which no two indices are the same, for example, 

x x , x 2 , ... x r) y r + x , ...,y n (a.) 

in terms of the remaining n variables, 

y^y^ '••yn x r+i> •••»#»» (p.) 

then, taking Q= — (X)+ {x x )y x + (x,)!?,+ . r + (a:,)y, 

(in which the brackets indicate that the variables of the set (a.) are to be expressed 

in terms of those of the set (£5.), so that Q is a function of the latter set), we shall have 

x^=# s trom ^=0 to i=r, 
^=-y,from./=r+l to,;=», 
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and J£+d£—°> as before * i 

but the equations corresponding to (10.) will not subsist unless X be homogeneous 
with respect to the r variables x 1} x 2 , ... x r . 

5. Let us now suppose that the function X contains, explicitly, besides the n vari- 
ables x x ,x 2 , ... x n , another variable t, and also rc constants a 1} a 2 , ... a n ; and that these 
last are contained in such a way that the n equations 

dX . dX dX 

dT= b »dT= b »-dT= b » C 11 -) 

would be algebraically sufficient to determine a 1} a 2 , ... a n in terms of b u b 2 , &c, x u &c. 
Then taking X b =-(X)+(a 1 )b 1 +(a 2 )b 2 + ..+{a n )b n 

(the brackets indicating that a 1} a 2 , &c. are to be expressed as above supposed), we 
shall have, by the theorems of arts. 2 and 3, 

dXi_ dX b _ dX„ 
db 1 — a "db 2 — a2 '-—db n - ani ( U -> 

and also, for all values of i, 

dXj_ dX__ , , 

dan~ ~dx— ~^ ; {16.) 

to which we may add 

dX ± __dX_dJ 
dt ~ dt~~ dt "*'> 

Now assuming the 2n equations (5.) and (11.), namely (for all values of i), 

dX_ dX 
dx—y* dai~~ bh 

we may suppose each of the 2n variables x 1} x 2 , ... y x , y 2 , .... to be expressed by means 
of them as a function of the 2n constants a„ &c, b lt &c, and t ; or, conversely, each 
of the 2ra constants to be expressed as a function of the variables x 1} hc.,y x , &c, and t. 
On the former hypothesis each of the variables x l} ...y u .... is given as an explicit, 
and on the latter as an implicit function of the single variable t, which we will con- 
sider as independent ; and total differentiation with respect to t will throughout this 
paper be denoted by accents, which will be used for no other purpose. Thus, p being 
any function of all the variables, we shall have 

dp , * ,, ,4 , , 

For the rest, we may, when necessary, distinguish the meanings of the various partial 

* Although these theorems, as stated in the text, are more general in form than those of the preceding 
article, they may, under another point of view, be considered as particular cases of them, and may in this way 
be best established. 
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differential coefficients employed, by referring to the hypotheses on which they are 
taken, and which I shall denote as follows : — 

Hyp. I. — The 2n variables a?„ x 2 , ...y 1} y 2 , ... expressed as functions of a l} « 2 , ... b lt 
b 2 , ... and t. 

Hyp. II. — The 2n constants a u a 3 , ... b 1} b 2 , ... expressed as functions of a?„ x 2 , ... y 13 
y 2 , ... and t. 

Hyp. III. — The n variables y 1} y 2 , ... y n expressed as functions of the n variables 
x„ x a , ... x n , the n constants a u a 2 , ... a n , and t (as by equations (5.)). 

Hyp. IV. — The n constants b 1} b 2 , ...b n expressed as functions of the n variables 
x 1} ... x n , the n constants a„ ... a„, and t (as by equations (11.)). 

6. Differentiating totally the equation (11.), 

dX_ 
da~ bi 

with respect to t, we obtain (observing that J^=^. by virtue of the conditions (5.)), 

da^da* l ^dai 2 ^ — + dot n ~ U 

( whei"e j^j &c. are taken on Hyp. III., art. 5. ). 
Now let (Z) be a function of x 1} ... x n , t, «„ ... a n , denned by the equation 

(A)= ~^> • ........... ('&•) 

the above equation then becomes 

fl- rf l! r ' + k'+ 4-^V 

If this equation be multiplied by x» and the result on each side summed with respect 

to i, it will be seen that the coefficients of x\, x 2 , &c. on the second side all vanish 
except that of x'j, which reduces itself to 1 (see art. 1, equations (1.) (2.)); so that 

we have 

d(Z) da x , d(Z) da^ ,d(Z) da* , 

da x dyj'da.% dyf* "'* da n dt/j ?' 

Now the expression on the left of this equation is equivalent to 

dZ 

d Vi 

if by Z (without brackets) we denote the result of substituting for a x , a 2 , ... a n in (Z), 
their values in terms of all the variables {Hyp. II.), so that Z is a function of the 
variables only. We have then, finally (writing i instead of/), 

*;=?,• ceo 

Again, we have (Hyp. III.) ,,:=§+**;+*<+.... 



dt ' dx x l 'dx. 



2 
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which, by (5.), (15.) and (16.), becomes 

, d[Z)_dZ d^ <zz % 

"'■ cfoj dyj dx-^dy^ dx^'"' ' 

. . . . , d(Z) dZ dZ dy l dZ dy % 

" dxi dxi dy x dxi'dy^ dxi "' 

(since (Z) would be derived from Z by substituting in the latter the expressions for 
y„ y 2 , ... Hyp. III.). And since -^=.-p-, &c, comparing the two equations last written, 
we obtain 

^=-1 «•> 

The system of 2rc equations (16.) and (17.) express the result of eliminating the 
2n constants from the equations (5.) and (11.) and their differential coefficients with 
respect to t. In other words, (16.) and (17.) are a system of 2« simultaneous differ- 
ential equations of the first order, of which (5.) and (11.), or again, the equations 
supposed in Hyp. I. or II., art. 5, are the 2n integrals. 

7. There are other remarkable relations between the partial differential coefficients 
of the expressions supposed in Hyp. I. and II., art. 5. For if we differentiate the 

equation -jz=b { with respect to a s {Hyp. I.), we obtain 

<? 2 X d*X dxj jFX_ dw^ 
daidaj'da i dx l daj'da i dx i dap""' ' ^ '' 

dli. dX 

which gives, putting &,■ for -t— and y t for -y— , 

(Kit &3Ci 

ddi ddi ddjdcii daj''"'dai daj *• '' 

(where j|> refers to Hyp. IV., J- 1 , &c. to Hyp. III., and J?, &c. to Hyp. I.V 

If then this equation be multiplied by -r~ {Hyp. II.) and the result summed with 

db- 
respect to i, the sum of the first terms is jJ- {Hyp. II.), and for the rest, the coeffi- 

dcci 

cient of j- reduces itself to unity, whilst those of the remaining terms vanish (art. 1 , 

equ. (1.), (2.)). Thus we have 

d h.— ^? na\ 

dy k —da i V s 'l 

(where the first side refers to Hyp. II., and the second to Hyp. I.). 

Now if we treat the equations 

dY_ dY_ 
dy— Xi > da—~~ bi 
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(see equations (6.) and (9.), putting a t for p in the latter) exactly in the same way, 
it is plain that the result may be deduced from (18.) by interchanging x and y, and 
changing the sign of b ; thus 

dbj dyu 
dxjt da/ 

Lastly, from the equations 

^.^-^ (see (12.) and (13.)), 

we should find in a similar manner 

da } - dxj, . 

dyit dbj 

and from the analogous equations (the existence of which is obvious) 

d yi -~ Xi > db t — Qi 
we should obtain £=-^' 

Collecting these results, and changing the indices, we have the system 

.......... (19.) 



dxi dbj dwi daj 

daj dt/i dbj dy. 



dyi dbj dyi daj 

daj dxi dbj dx t 

in each of which equations the first member refers to Hyp. I., and the second to 

Hyp. II. (art. 5.) ; and it is to be remembered that there is no relation between the 

indices of the variables and those of the constants, so that the case of i=j has no 

peculiarity*. 

8. Let $, A be symbols denoting two distinct sets of arbitrary and independent 

variations attributed to the 2n constants ; then the equations 

dX dX. , 

dxr^y* lai~~ bi 

give SX=%>;,.+«^i) ; 

and if the operation A be performed on each side, we have 

ASX=2 i (A 3 /M+A« i ^ i ) 
+%(yA^ i +aAto i ). 

* It is remarkable that each of the equations (19.) is also true on a different and separate hypothesis, as is 
apparent on inspection of the four different sets of equations, 

dX_ dX b _ dY_ dY b _ 



dX. 7 dXh dY » dYi 

-y-=Oi, -—-=ai, —~=—bi, — — = — a, 
aui dbi ddi dbi 



(see the preceding articles). 



DIFFERENTIAL EQUATIONS OF DYNAMICS, ETC. 81 

If from this we subtract the corresponding equation obtained by inverting the order 
of the operations A, d, remembering that &bu=l/\u, we obtain 

XfaAyi-tefyd+liQaAb-bafa)^* • ( 20 -). 

(The use here made of the double operation AS, is clue in principle to Mr. Boole. 
See his demonstration of a well-known theorem of Lagrange, of which the equation 
(20.) is a more general form-|~). 

If in this equation we suppose $x i} ly i} &c. to be expressed in terms of la { , hfy, &c. 
{Hyp. I.), and A^, A& i5 &c. in terms of A# i} Ay if &c. (Hyp. II.), and compare the terms 
on the two sides, it is easy to derive the relations (19.). I preferred however to 
deduce them by a more direct method. 

9. If Xf be expressed in terms of the 2ra constants and t (Hyp. I.), and then each 
constant be expressed in terms of the variables (Hyp. II.), the result is an identical 
equation. Differentiating then with respect to x i} x J} y k , we obtain the three equations 

dxi da x d%i db x dxi da% , dxi db 2 
«, d%i'db 1 dxi da^ 



UX( ««j WDi UU X UX{ UUq UXi «W 2 „ 

da, dxi'db x dxida% dxi'db^ dx,' 



dx, da x dxi db x dxi da, 2 dx f db% „ 

da x dxj'db x dxj'da^ dxj'db^ dxf* 

dx { da x dxi db x dx t da% dx t db^ . ^ 

da x dy k ~*db x dy^'da^ dy k 'db q dy k 

Three similar equations may be obtained by treating y i in the same way. And if we 
apply to these six equations the transformations given by the system (19.), art. 7, 
the resulting theorems may be comprehended in the following statement. 
If p, q be any two of the 2« variables a?„ ... x n , y x , ... y n , then 

„/dpdq_ dp_^\_^fdb i da i dhdaA 

Z, \db i da i "~da i db { )— 2 '\dpdq dq dp J — ^ l > or ~ u ' • • • K**-) 

according as p and q are or are not a conjugate pair, i. e. a pair of the form x p y^. 
(The value +1 belongs to the case in which p=x J} q=y s , and —I to the converse.) 

Here^j and q are a determinate pair of variables, and the summation refers to the 
constants, extending to the n conjugate pairs. 

More important however are the converse theorems obtained in a perfectly similar 
way by expressing a t , or b ( in terms of the variables (Hyp. II.), and supposing the 
variables to be again expressed in terms of the constants and t (Hyp. I.). Differenti- 
ating the resulting identical equation with respect to a { , %, b { , b s , and applying the 
transformations (19.) as before, we have, putting h, k for a determinate pair of con- 
stants, 

/dh_ dk_ dh_ dk\ _~ (dy±d?h Ax% AyC\ ,, _ n (i)0 \ 

2 '\dy i 'dx~dx i dy i ) — z '\dhdk~dhdk) — '^ l ' 0l '— K} ' ' ' ' v"'' 

* This might be written 

S ; 5(a;,-,y,) + 2^(a,-, b,)=Q. 
See the notation proposed in art. 1 . 

f Cambridge Mathematical Journal, vol. ii. p. 100. 
MDCCCLIV. M 
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according' as h, k are or are not a conjugate pair, i. e. of the form a j} bj. (The value 

+ 1 belongs to h=a p k=b p and —1 to the converse.) 

According to the notation proposed at the beginning of this paper, the above 

formula may be written 

d{h,k) _^ d{y { , Wj) , _ 

Z 'd{ Vi , w t )~^ d(h, k) ~^ J ' 01 - U - 

By a usual and convenient abbreviation, the sum 

s «*(*» *) 

may be denoted by the symbol* [h, Jf\, We have then, by (22.), 

[a,,i,J = — [6 a,] = l [«,-, i y ] = [a { , aj = [>, ^]=0, (23.) 

j being different from i ; and, obviously, 

[a,-, a,] = [i &J=0. 

Now let/, g be any two functions whatever of the 2n constants a„ &c. b 1} &c. ; when 
the latter are expressed in terms of the variables {Hyp. \l.),f, g become also functions 
of the variables ; and if h, k represent, as above, any pair whatever of «„ &c, b 1} &c, 
we have (see art. 1.) 

d(f,ff) = Mf,ff) d(h,k) ] 

<%« %i) \d{h, k) ' d(y { , os t ) J ' 

the summation referring to h, k, and extending to every binary combination. 
If, now, we sum each side of this equation with respect to i, we obtain 

[/, ri =s{[M].«!} (24.) 

(the summation referring as before to h, k). But, by (23.), [h, k~] is unless h, k be 
a conjugate pair, and then it is +1 ; so that (24.) becomes simply 

U,S\^^y (25.) 

an equation which, written at length in the common notation, is 

1 \dyi d%i dw ( dy { J \dai dbi db t da { J' 

The expression on the right being a function of the constants a l} &c, b x , &c. only, the 
equation (25.) expresses obviously the following theorem. 

If f=<P{ x i> x » — x n, y» y» ••• y« t) 

g=Tp(x l} x„ ... x n , y 1} y„ ... y n , t) 
be any two integrals of the system of simultaneous equations (16.), (\7-), art. 6, then 

* Poisson employs the notation {h, k), which would have led to confusion if adopted here. Lagrange (in 
the Mic. Anal.) uses [h, k], but with a different signification. See below, note to art. 34. 
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the expression [/, g], or 






d<p dty 
d%i dyA' 



is constant ; i. e. it becomes a function of the arbitrary constants only, if for x t , &c, 
y i} &c. be substituted their values in terms of the constants and t. 

In the case in which (16.) and (17.) represent the dynamical equations, this is 
identical with the remarkable theorem discovered by Poisson. We shall have occa- 
sion to return to it presently. 

10. If we treat the equations (21.) of the last article exactly in the same way as 
we have treated (22.), putting w, v for any two functions whatever of the 2n variables 

x 1} x%, . . . x n , y 1} y 2 , . . . y n , 

/. i ^ d(u, v) 

and comparing this with the theorem (25.) of the last article, we see that both may 
be included in the following general enunciation : — 

If u, v be either (1) any two functions whatever of the 2ra constants a lt &c, b i} &c, 
or (2) any two functions whatever of the 2» variables x 1} &c, y u &c. (not containing 
t explicitly), then 

du do 



, d(u, v) ^ 
H d(xi, y { )' 



or 



-,\du dv du do * „„ „„, „.. 

~' i \dyi dxi dwi dyi'dbi da ( dai db { \ ' 

d(u, v) d(u, v) 
dfyi, %i)' d{b, 



m=o. 



(26.) 



(When u, v represent functions of the constants, the differential coefficients in the 
first term are taken on Hyp. II. ; and, when functions of the variables, those in the 
second term on Hyp. I. (art. 5.)). 

This property depends, as will be seen, solely on the relations (5.), (11.), arts. 2, 5, 
which are the only assumptions that have been made in deducing all the preceding 
propositions. 

11. There are similar theorems in which the summation refers to the numerators of 
the differential coefficients ; but as these are less remarkable, and moreover are 
deducible immediately from the equation (20.), art. 8, 1 shall omit them. 

12. Theorem. — I proceed now to establish a theorem which may be considered as 
the converse of that expressed by (23.), art. 9. 

Let x l} x 2 , ... x n , y u y 2 , ... y n be 2n variables, concerning which no supposition what- 
ever is made, except that they are connected by n equations 



"l — < Pl\ X l> *25 ••• ^nJ Vli y%> ■•• Vn) 

a !i =-<p2{x 1 , x 2i ... x n , y x , y 2 , ... y n ) 



(a.) 



a,, a„ ... a n being n constants. The functions on the right may involve explicitly any 

m2 



*1> "23 
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other quantities whatever, except a„ &c. It is assumed that these equations are 
algebraically sufficient to determine each of the n variables y 1} ...y n , as a function of 
the other n variables x u ... x n and the constants. Then the theorem in question is 
as follows : — 

If, by means of the equations (a.), the n variables y„ ... y n be expressed as 
functions of a?„ &c, then in order that the conditions 

dy 1 _dy l 

Ujvb a (XiXi j 

may subsist identically, it is necessary and sufficient that the expression [a t) a,] 
(denned as in art. 9.) shall vanish for every binary combination of the n equations. 
This may be proved as follows : — 

Putting- h, k for any two of the constants a 1} a 2 , &c, let h=<p(x l , &c, y 1} &c.) repre- 
sent one of the equations (a.) above written. If in this equation the values of «/„ ...y n 
be expressed, as above supposed, in terms of #„ &c, a„ &c, it becomes identical 
Differentiating it, on this hypothesis, with respect to #,., we obtain 

JA <tt <fe/, <?A rfy 2 <tt dy n 

dx i 'dy l dxi'dy^ dxi''"'dy n dx, ' 

and in like manner 

M.^dyidk_dy2 dk dy n _ 

d%i ' dy x dx{ ' dy 3 dxt''"' dy n dx ( ' 

dk dh 

and if we multiply the first of these equations by j- and the second by j- and sub- 
tract, there results an equation which may be written as follows : — 



dh dk dh dk (dyj/dh dk dh dk\~\ 

dyi dxi dxi dy t AdxAdyj dyi dy ( dyjj J ' 



or, putting now a p) a q instead of h, k, and employing the same notation as before, 

d{a p , a q ) __^\dyj d{a p , a g ) \ 
d(yi,Xi) \dxi d{yj, yi)y 

If now the terms on each side be summed with respect to i, the result on the first 
side is [a p) aj ; and observing that on the second side the term multiplied by -P will 

only differ in sign from that multiplied by-^» we shall have 

K,„,]=^{(| : -f;)*^} m 

the summation on the right extending to all binary combinations i, j. Suppose this 
equation to be written at length, and then after multiplying each side by 

<%-, y,) } 
d(a p , a q ) 

let the sum be taken with respect to all the binary combinations p, q. It follows 
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from the theorems of art. 1, that the coefficient of 

dy 1 _dy ± 

dtXfg UuXf 

on the right will reduce itself to unity, and that of each of the remaining terms to 
zero ; so that we shall have, writing nowj, i for r, s, 

i-ir^h^-^ (28 -> 

In order then that the expression gz~d snot, ld vanish identically for every binary 

combination of indices, it follows from (28.) that it is sufficient, and from (27.) that it 

is necessary, that each of the — ~ — - terms \_a p , a q ~\ should vanish,, and vice versd. It 

will be observed that the terms [a p , aj cannot vanish otherwise than identically, 
since they do not contain any of the constants a„ a 2 , &c, and it is by hypothesis im- 
possible to eliminate all these constants from the equations (a.). It follows then that 
when the conditions [a p , aj =0 subsist, the values of y u ... y„ expressed as above, are 
identically the partial differential coefficients of a function of x 1} ... x„, a u ... a n . 
We have thus established the theorem enunciated at the beginning of this article. 
13. The preceding theorem may be made somewhat more general as follows : — 
If we divide the 2n variables into any two sets of n each, so that no two in the 
same set are conjugate (as for instance 

X l} X 2 , ... X f , ffr+li ••• Vn 

and denote one set by 
and the other by 



SWJ 



SlJ §2? • • • Sra 

taking the + or — sign according as *? £ represents y t or x ( , it is obvious that the ex- 
pression 2 4 4 p ' j! is identical with [a p , «J ; and therefore whenever all the terms 

[a p , «J vanish, if the set r\ x , %, ... q n can be expressed by means of the equations (a.) of 
the last article, in terms of ?„ §,, ... l n , a„ a a , ... a„, their values will be the partial 
differential coefficients with respect to £ 1} | 2 , ... |», of a function of these variables and 
of the constants. 

14. Theorem. — If of the system of 2n simultaneous differential equations of the first 
order 

Xj -j — 5 ...., x n — — 



(I.) 

(where Z denotes any function of x u ... x n , y 1} ... y n and t, and accents denote as usual 



>__dZ •—_ d ?L 
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total differentiation with respect to t) there be given n integrals, involving n arbitrary 
constants a„ ... a n , as 

the remaining integrals may be found, whenever the g conditions [a i} aj\ =0 are 

satisfied. 

For let y u y 2 , ... y n be expressed, by means of the given integrals, in terms of 

x x , ... x n , a 1} ... a ni t. 
Their values so expressed will satisfy (art. 12.) the conditions 

Let (Z) represent the result of substituting in Z these values of y x , ...,y n , so that 
(Z) is a given function of*,, ... x n , a u ... a n , t. We shall have 

d(Z)_dZ dZ dyi dZ dy^ 
dxi dx~ r dy x dxi^dy^ dxi'"' 

which the equations (I.) and (&.) reduce to 

dx t — y^dxp ^dx^ -" 
but «'-* + ^±^'4. 

consequently l^^-ft (*) 

Looking now at the assemblage of equations (b.), (c), we see that they express the 
following proposition :■ — 
The values of y u y„ ...., y., -(Z), 

are the partial differential coefficients with respect to x 1} x 2 , ... x H , t, of one and the 
same function. Let this function be called X ; we have then 

dX ^ X /rr\ /it ^ 

^■=^^-=-( z ); (»•) 

and since y 1} ...,y n , (Z) are given functions of x u &c, a x , &c, t, the function X can 
be found by simple integration. 

Let us then suppose X to be known, and let us take the total differential coefficient 

with respect to t, of -^ ; we shall have 

(dX\'_d^X_ , d*X , , i 2 X #a? , , 

\dai) dciidt ddidx x ' da,idx% 2 
which, by virtue of (I.) and (II.), becomes 
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eki dy t da { ~^ dy a dai* " " 

(since (Z) is derived from Z by introducing the values of y l} ...y n , in terms of x 1} &c, 
a„ ... a„), hence the second member of the preceding equation vanishes, and we have 

(£)'=»> 

so that -T- is constant, and we may write 

S= 4 » • • <»'•> 

and b t is an independent arbitrary constant, as it is easy to prove ; it is however 
unnecessary to do so here, because we have in fact already proved it in showing that 
the elimination of a„ ... a„, b 1} ... b„, from the system of equations (II.), (III.), leads 
to the differential equations (I.) (see art. 6.). The n equations (III.) give therefore 
the remaining n integrals of the system (I.), of which (II.) and (III.) together are the 
complete solution. 

The system of equations (II.), (HI.) being the same as that discussed in the pre- 
ceding articles, all the conclusions there obtained will continue to subsist. 

15. Suppose the expression for Z (see the last article) in terms of the variables is 

£t^~/(Xi, x i} ... x„,y if y S) ...y„, t), 

Z is changed into (Z) by the substitution of jg- f° r Vu & c * » an d s i nce ~^r is (identi- 
cally) = — (Z), the equation 

f+y(x„,„...x„g,...,g, ( )=o (X., 

is a partial differential equation satisfied by the function X. 

We have thus arrived, by an inverse route, at the point from which Sir W. Hamil- 
ton's theory, as improved by Jacobi, sets out. 

Jacobi, namely, has shown (by a demonstration immediately applying only to a 
particular form of the equation (X.), but easily extended), that if X be any "com- 
plete" solution of the equation (X.), that is, a solution involving (besides the constant 
which may be merely added to X) n arbitrary constants a lt o 2 , ... a n , in such a way 
that they cannot be all eliminated from the ra+1 equations obtained by differentiating 
X with respect to x u ... x n , t, without employing all those equations, then X possesses 
the properties of Sir W. Hamilton's " Principal Function," or in other words, gives 
all the integrals of the system (I.) by means of the system (II.), (III.). It will be 
desirable briefly to indicate the mode in which this demonstration may be made to 
apply to the general form (X.). 

Assuming that a complete solution X, of that equation, is given, put T-=y <; then 
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differentiating the equation (X.) with respect to x { , and employing the equations 

G/Xq (IX p 

dt'dx i 'dy 1 dx^'dy^ dx^'' ' 

on the other hand, taking the differential coefficient of y t with respect to t, without 
assuming anything as to the nature of the relations between t and the other variables, 
we find 

#«— dt + <fe 1 a?1 " t \fe/»T"" 

and adding to this the preceding equation, 

"'' dxi dx^ i dy l J'^dx^\ 2 dx % J'"" 
from which it follows that the n assumptions 

x *-d Vi 
would involve the n further equations 

,__df 
Vi— dssi 
Again, the n assumptions 

dX . 

would give, by combining the n equations obtained by differentiating totally with 
respect to t, viz. 

ddidt ' da i dx l 1 ' da,dx^ 2 >"" ' 
with the n others obtained by differentiating the equation (X.) with respect to a { , viz- 

d*X df d*X df d*X 

daidt dy 1 da^Xi dy% datdx^'"' ' 

the n following, namely, 

d?X / , Jf\ _^X / , _f \ 
da i dx 1 \ Xl ~dyJ + da i dxX X2 ~dyJ + °' 

from which it follows either that a?J= -/-> or that the determinant formed with the « 2 

expressions t—t-j or — (-r-), vanishes; but this last condition would express, as is 

well known, the possibility of eliminating the n constants a„ a 2 , ... a n from the n 
equations 

which would contradict the assumption that X is a complete solution of the equa- 
tion (X.). 
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Finally, then, if X be a complete solution, the assumptions -r— =ij involve as a 

consequence the relations x'^-J--, and these again involve y\=.—J~, where y t stands 

. dX 
tor j — 

axi 

In thus applying Jacobi's demonstration I have slightly altered its form, in order 
to bring more prominently into view the necessity for X being a complete solution. 

16. It is obvious, from the considerations given in art. 13, that instead of the 
equation (X.) of the last article, we might employ any one of the analogous equations 
obtained by distributing the variables as explained in the article referred to, and then 

writing ^f for % in the expression for Z. The function Q will be a " principal func- 
tion." In particular, if we take the equation 



dY JdY dY /I n 



any complete solution will give the integrals of the differential equations (I.) by means 

of the system 

dY_ £Y, 

dyi " da.i~~ *' 

The whole number of partial differential equations from each of which a " principal 
function" can be obtained, will obviously be 2 n . The relations between these dif- 
ferent principal functions will be apparent from the conclusions of art. 4*. 

17- If #u ^ 2 > ••• x n represent all the independent coordinates (of whatever kind) in 
any ordinary dynamical problem, and T the expression for the vis viva-f- in terms of 
x l} &c, x\, &c, the equations of motion are, as is well known, 

ydx'i/ dx t dxi' ■■■•••■•• • \ ■) 

where U is a function of x 1} ... x n , which may also contain t explicitly, but not x\, &c. 
Lagrange, to whom these formulae are due, was also the first to employ the expressions 

jrp 

—i as new variables, instead of*-. But Sir W. Hamilton first showed that this sub- 

dxi 

stitution (putting —,=y i ) would reduce the n equations (T.) to the 2n equations of 

the first order of the form (I.), art. 14. His demonstration, however^, depends upon 
the circumstance that T is, in dynamical problems, necessarily homogeneous with 
respect to x\, .... x' n , and I am not aware that any other case has hitherto been con- 
templated. 
The investigations of the preceding articles will however enable us to apply a 

* Compare Sir W. Hamilton's expressions, Philosophical Transactions, 1835, p. 99, art. 5. 

f I here adopt, what I hope will be universally adopted, the suggestion of Coeiolis and Professor 
Hklmholtz, that the definition of vis viva should be half the sum of products of masses by squares of velocities. 

X Philosophical Transactions, 1835, p. 97, art. 3. 
MDCCCLIV. N 
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similar transformation to the equations (T.), in the case in which no limitation is im- 
posed upon the form of the function T, as I shall now proceed to show. 

18. Putting T+U=W, we shall have (since U does not contain x\, &c.) 

m=z ■ (wo 

dW , 

Let — T=y { ; then if we take 

z=-(W)+(a?;) yi +(4)y,+ ...+(^)y. • (V.) 

(where, in the terms enclosed in brackets, x\, x\, &c. are to be expressed in terms of 
y y , y 2 , &c, x 1} x 2 , &c), we shall have, by the theorems of the former articles (see equa- 
tions (6.), (8.), (9.) of arts. 2 and 3, putting- x- instead of x ( and x t instead of p, in those 
equations), 

(«•) 

and f" 

i 

so that the equation (W.) becomes 

03.) 

and (a.), (/3.) are of the form in question* ((I.), art. 14.). Thus, so far as the appli- 
cation of any methods of integration, founded upon the preceding principles, and the 
theories of Sir W. Hamilton and Jacobi, to the system (T.), art. 16, is concerned, 
there is no restriction to the form of the function T. This extension is probably at 
present of no practical importance, but may perhaps be thought of some interest in a 
purely analytical point of view. 

19. Returning now to the suppositions and conclusions of art. 14, let us further 
suppose that Z does not contain t explicitly, so that 

by virtue of the system (I.) ; in this case 

Z=A (h.) 

is one of the integrals of the system, and if we suppose this to be one of the n given 
integrals from which the principal function X is to be found, so that 

n, a l) # 2 , ... ct n _ 1 
are now the n arbitrary constants, and the conditions 

[ff t -, %]=0, [h, oj=0 
subsist, it is plain that we shall have 

(Z)=A, 

* In the case in which T is homogeneous and of the second degree, in x\, x' 2> ... x n , it is obvious that the 
expression for Z reduces itself to 2(T)— W, or (T) — U. 
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since the expression for Z must reduce itself identically to h when the values of y^ ...y n 
obtained from the integrals are substituted in it. Hence 

dX 



and therefore X=— ht-\-Y, 

V being a function not containing i 

to be fouhd from the n expressions 



V being a function not containing t explicitly. We have then — =— - , so that V is 



dj__ 
dx~ Vi ' 
Lastly, the n remaining integrals will be 

dh ' da,i * 

(r representing the arbitrary constant conjugate to h) ; and, substituting in these the 
above expression for X, we obtain 

dV ,, dV , , nn . 

M= t+r > i^r h (29 - } 

The function V now satisfies, and may be defined by, the partial differential equation 

/(*»...*. s ->...» a j=A, ( v.) 

where/(«„ ... x n ,y n ...y n ) is the expression for Z in terms of the variables. 

This, in dynamical problems, is the case in which the so-called " principle of vis 
viva" subsists. I shall, in the rest of this paper, use h exclusively in the above signi- 
fication, and call it, whether actually referring to a dynamical problem or not, the 
"constant of vis viva," whilst the integral Z=h may be called the "integral of vis 
viva." 

20. When the 2n integrals of the system of differential equations (I.), art. 14, are 
expressed in the manner which has been explained, it follows from the conclusions of 
former articles, that when these integrals are put in the form 

o«= &(#» .-., x n , y 1} ...y n , t) 

i 4 =^,(ar„ ..., x„y lt ... y n , i), 
the conditions [<*„ J,] = 1, [a„ ft,] =0, [b i} bj] =0 will subsist, as well as [a l} «,] =0. I 
shall call any system of 2n integrals in which these conditions are fulfilled, a " normal 
solution," or a system of "normal integrals," whilst the 2ra arbitrary constants con- 
tained in such a system may be called "normal elements." Any pair a t , b i} may be 
called (as before) conjugate elements. In the case considered in art. 19, h and r are 
conjugate elements, these letters being used instead of a, b, merely from obvious 
motives of convenience. 

It has been one principal object of these investigations to ascertain what advantages 
could be gained— either for the actual integration of a system of equations of the 

N 2 
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form (I.), or for the transformation of known solutions into forms convenient for the 
application of the method of variation of elements — by making the discovery of prin- 
cipal functions depend upon that of n integrals satisfying given conditions, rather 
than upon the solution of a partial differential equation. Having now prepared the 
way for this inquiry, I shall proceed with it in the following section. 

Section II. 

21. Theorem.— If p, q, r be any three functions whatever of the 2ra variables 
%i, •-, x M y» •'•?/«> then 

[[?» ?].»•]+[[?» r ]»y]+[[r,p],?]=0 (30.) 

(The symbols have the same signification as in the last section. See art. 9.) 

This may be proved as follows. It is evident that if the above expression were 

developed, each term would consist of a second differential coefficient of one of the 

functions/), q, r, multiplied by a. first differential coefficient of each of the other two. 

Consider then the terms in which p is twice differentiated ; these will be of the 

three forms 

d*p dq dr d?p dq dr . d 2 p dq dr 



?-, and-; 



dxidyj dxj dyt dxidxj dyt dyj dy t dyj dx ( dxj 

each of which will arise from the first and third terms of (30.) only. (It is to be ob- 
served that i may ==/'.) 

Now if we examine each of these forms, we see easily that for every term arising 
from the first term of (30.), there is a similar term with the opposite sign arising from 
the third term of (30.) ; and since a similar proposition would be true of the terms in 
which q, r, respectively, are twice differentiated, the whole expression on the left of 
the equation (30.) vanishes identically. The theorem is therefore established. 

It is obvious that p, q, r may contain, explicitly, any other quantities (as t) besides 
the 2n variables with respect to which the differentiations are performed. 

Let £ represent, either, one of the 2rc variables *„ &c, t/ M &c, or any other quantity 
whatever, explicitly contained in p and q. It is evident that we shall have 

*[» *>GH+ [»|] (»>•) 

22. Resuming now the consideration of the 2ra simultaneous differential equations 
discussed in the first section, namely, 

Xi -ij: *-- s? ( L > 

we shall be enabled, by means of the theorems (30.), (31.) of the last article, to give 
a very simple and direct proof of the proposition indirectly demonstrated in art. 9. 
For let u be any function whatever of the variables x„ &c, y„ &c, t ; then 



*=§+K^+t4 
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and if the values of x' t , y\ given by (I.) be substituted in this expression, it becomes 

rf=5 + [Z f «]. (32.) 

Let w=[/>, q], then (making use of (31.)) 

[?. ?]'= [|> f]+ [r, §]+[z. [*.»]]• 

Now suppose that, by virtue of the differential equations (I.), the values of p and g 
are constant ; or, in other words, that 

p=<p(x 1} &c, y„ &c, *) 

g , ='4'(*i 5 &c.,y I , &c, 

are any two integrals whatever of the system (I.) ; p, q representing two arbitrary 
constants. The equation p'=-0 gives (see (32.)) 

t=-c z .rf. 

hence [§, }] = -[[Z, P ], «]=[,, [Z, ,]], 

In like manner fp, ^ |=fp> [? 3 Z]T 

Thus the expression given above for [jp, §•]' becomes 

[ P , q y=[p, [,, z]]+[^ [z, p]]+[z, o, 9 ]], 

which is identically equal to 0, by the theorem (30.). Consequently, for any two 
integrals p and q, 

[p, q]= constant (33.) 

This theorem, as has been already mentioned, was discovered, in the case of the 
dynamical equations, by Poisson ; and the fact that he was able to arrive at it 
through so long and complex a process as that which he has given in his first memoir 
on the Variation of Arbitrary Constants*, must be looked upon as a remarkable 
instance of his analytical skill. I am not acquainted with any attempt to simplify 
the demonstration, except that of Sir W. Hamilton -f~ ; in fact it is probable that no 
material simplification was attainable without the help of the transformation of the 
differential equations to the form (I.), towards which Poisson (as Jacobi has remarked) 
only made a first step. Sir W. Hamilton's demonstration may certainly be con- 
sidered simple as compared with that of Poisson. That which I have given above 
will, I hope, be regarded as a further improvement. 

23. In what follows I shall use such expressions as "the integral c," as an abbre- 
viation for "the equation c—^x^ &c, #„ &c, t)." 

* Journ. de l'^lcole Polytechnique, torn. viii. t Philosophical Transactions, 1835, p. 108-9. 
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It is of course understood that the function on the right contains neither c nor any 

other arbitrary constant explicitly. 

Let then /, g be any two given integrals of the system (I.). It has been shown 

that we shall always have 

[f, g] = constant (K.) 

But this equation may be true either (1) identically, or (2) not identically. In the 
first case the expression [/, g~] may either be identically =0, or it may reduce itself 
identically to a determinate constant, which might always be made unity by multiply- 
ing one of the integrals by a factor. (In the case of a " normal system" of integrals 
(art. 20.), it has been seen that every binary combination gives either or 1.) But 
if the above equation (K.) be not identically true, so that [/, g~] obtains a constant 
value only by virtue of the differential equations, then the constant on the right of 
(K.) is an arbitrary constant, and that equation is itself an integral. But here again 
there are two cases ; for the function [/, g~] may be only a combination of the func- 
tions on the right of the two integrals f } g ; and then (K.) is not a new integral, but 
only a combination of the two given ones ; or, on the other hand, [/, g] may be a 
function independent of,/', g; and then (K.) is really a new integral, which cannot 
be produced by merely combining the other two. Thus it appears that Poisson's 
theorem may in some cases lead to the discovery of new integrals, when two are known. 
On this subject, and others connected with it, I refer to the interesting memoir of 
M. Bertrand in Liouville's Journal (1852), "Sur Fintegration des equations diffe>- 
entielles de la Mecanique." 

24. Let c x , c 2 , ... c m be any m integrals, and let f, g be any two functions of the m 
constants c„ c 2 , ..., c m , so that/, g are also two integrals; and considering /, g as 
functions of c„ ..., c m , and, through them, of the variables, we have exactly as in 
art. 9, equation (24.)_, 

w*H$£iM (L -> 

the summation extending to all binary combinations of the m constants c„ &c. If 
then we suppose k lt k 3 , ..., k m to be m functions (such as h, k) of the m constants 
c„ ..., c m , we shall have for any pair k p , k g , 

ft^Htrtf^} (**■) 

(the summation referring as before to i,j) ; and the inverse equations (obtained either 
by considering c„ &c. as functions of fc„ &c, and reasoning in the same way, or by 

die c) 
multiplying the above equation by ^A " -p. and summing with respect to p, q) will be 

h%MfefrR.«j} < 35 -> 

(the summation referring to p, q). 
This inversion can only fail in the case in which the equations expressing #„ &c. in 
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terms of c 1} &c. are not all independent ; a supposition which we exclude, in order 
that Ar„ ... k m may represent m distinct integrals. 

The equations above written lead obviously to the following conclusions:— 

(1.) If/ be a given function of the m constants c 1? ... c m ; then the determination 
of another function g, such that [/, g~\ =0, depends in general upon the solution of 
a linear partial differential equation of the first order. 

(2.) It is impossible that the conditions [ft,, A:J=0 can exist for every binary com- 
bination of h u ..., k m , unless [c ( , c?]=0 for every binary combination of c l} ..., c m . 

25. As an illustration of the first of these conclusions, we may take a case which 
actually occurs in many dynamical problems. Let c„ c 2 , c s be three integrals, such 
that 

|_C 2 , C 3 J=C„ [_C 3 , CiJ— C 2 , [Ci, c 2 J=c 3 , . . . . . " . . (c.) 

and let it be required to find a function g of c x , c 2 , c 3 , such that [c„ jr]— 0. The 
equation (L.) of the last article gives, if we put/=c 1} and introduce the conditions (c), 

3 -z C^-; — V. 

<fc 2 

The solution of which is 



'<fc 2 dc 3 



g=^(<f a +cl), . . , (g.) 

■>p being an arbitrary function (which may evidently also contain c t in an arbitrary 
manner). 

If, instead off=c 1 , we put/=p(cHc2+c|!)> it will be found that the expression on 
the right of the equation (L.) vanishes identically; so that in this case, if g be any 
arbitrary function of c„ c 2 , c 3 , the condition [/, <?]=0 will be satisfied. 

26. If o„ a 2 , ... a„, b u b 2 , ... b K be a system of normal elements (art. 20.), we have 
(equation (25.), art. 9.) 

U '» J '<%, 6.) 

where /, g represent any two functions of the elements, or in other words, any two 
integrals whatever. If in the above equation we put successively f=a ( , f=b i , we 
obtain 

[»»J]=|. [*.»] = -£ (36.) 

In the case where the principle of vis viva subsists, we may suppose the constant of 
vis viva, h, to be one of the elements. In this case (see (29.), art. 19.) the element 
conjugate to h is r, and t appears in none of the integrals explicitly, except one, 
namely, the integral conjugate to A, which is 

If, then, g be any integral whatever, not containing t explicitly, it cannot contain r, 
since any combination of the normal integrals involving r, will involve it in the form 
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r+t. Consequently, for every such integral we shall have, by (36.), 

[g,h]=0, ........... (37.) 

since -f=0. 

ax 

This particular consequence of the formula (36.) follows also immediately from (32.), 
art. 22, since the equation g'=0 gives, by (32.), [Z, g}=®, and in this case Z=h, so 
that [Z, g] = [A, g]. In this manner the theorem expressed by (37.) has been already 
obtained by M. Bertrand. 

Examples of the preceding Methods. 

27. I shall now exemplify the principles which have been explained, by applying 
them to two of the most familiar as well as important problems of dynamics. First 
then let it be required to obtain in a normal form the integrals of the differential 
equations which determine the motion of a material point, acted on by a force ema- 
nating from a fixed centre and depending only on the distance. 

Taking the centre of force as the origin of a system of rectangular coordinates, let 
m be the mass, and x, y, z the coordinates of the moving point. Then 

T=lm(x>*+y<*+z% 

and U (see art. 17.) is a given function of r, say <p{r) t where r i =x*+y s +z i . Let us pot 

dT dT dT 

7i2= U > ^= y ' d?= W > 

so that, referring to the notation used in the preceding pages, we have 

x, y, z instead of x n x 2 , x 3 

u, v, w instead of y„ y iy y 3 . 
Moreover, u-=mx', v=my', w=mz'. 

Hence we obtain Z=(T)-U=w- 1 (« 2 +t> 2 +«> s )— <P(r), 

so that the integral of vis viva, or Z=h, becomes 

-(w 2 +t> 2 +tt> 2 )-p(r)=A ; 

and the three integrals which express the conservation of areas become 

yw — zt>=c, 

zu — xw—c 3 

xv—yu—c 3 . 
These integrals are immediately seen to satisfy the conditions 

from which it follows (see art. 25, the result of which is obviously unaffected by the 
negative signs), that if we take h— (cj+c5+cj)*, the condition [c 3 , ft] =0 will be satisfied 
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(as is easily found to be true) ; and since neither of the integrals c a , k contain t ex- 
plicitly, the conditions [h, c 3 ]=0, [h, k~]=0 will subsist also (art. 26.). Hence it 
follows that if we solved algebraically the three integrals h, c 3 , k so as to express u, 
v, w in terms oi x, y, z, their values would be the partial differential coefficients of a 
function V, from which the three remaining integrals could be found (arts. 12 and 19.). 
But it is more convenient to adopt a different system of coordinates. Reverting 
then to the primitive form of the three integrals which we have chosen, and writing c 
instead of c 3 , we have 

T-U=h (i.) 

m(xy' —yx')=c „ (ii.) 

w 2 (rV 2 +y 2 +3' 2 )-ry 2 )=& 2 . . (iii.) 

28. Let us now employ, instead of x, y, z, the three coordinates g, 6, z ; where z is 
the same as before, g is the projection of r on the plane of xy, and 6 is the angle 
between g and the positive axis of x. We shall thus have 

f-{-z*=r*, x=gcosd, y=g sin d, 
and T=lm( g '*+f0 2 +z»). 

¥ , dT dT dT 

Let d?= u > W= v > d7= w 

(where u and v have now a new signification), then 

,i « at— v „i— w . 



* *n 



:J Z —- 



m mg' m 

and the three integrals at the end of the last article become, after obvious reductions, 

±(*+$+*)=h+9(r). . (i.) 

v=c (ii.) 

(gw—zuy+- i v 2 =k i (iii.) 

The conditions [h, c]=0, [h, k~]=Q, [c, /c]=0 continue to subsist with reference to 
the new variables ; the two former necessarily, because (ii.) and (iii.) do not contain 
t (art. 26.), and the third actually, as is seen on trial {not accidentally, as will be 
shown hereafter). 

We know, therefore, that the values of u, v, w, found from these equations, will 
be the partial differential coefficients with respect to g, 6, z of a function V of these 

latter variables. 

<? 
The two first give u 2 -\-w 2 =2m(h+<p(r))— -5; 

and if we multiply this by f*-f 2 2 =r 3 , and subtract (iii.), we obtain (introducing the 

condition (ii.) 

(gu + zwf= 2mr !i (h-\-<p(r))- k\ 

MDCCCLIV. O 
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Lastly, if this be combined with (iii.), the following expressions are found for u and w : 



=i[2mr 2 (A+^(r))-^}*+^{^-^c^ 



(in which it is to be remembered that r 2 =z 2 -\-f) ) and if to these we join the equa- 
tion (ii.), the values of u, v, w are explicitly given in terms of the conjugate variables 
§, 6, z. We have then (art. 19.) 

\ =\(udg-\-vd&-\-wdz) ; 

or, substituting the above values, 

V=c^+J{^^(2mr 2 (A+^(r))-P)H ? ^^(^-->) i ]. 

The term under the integral sign is easily seen to be (as we know a priori it must be) 
a complete differential. It is convenient however to transform it thus. First, we 
have gdg-\-zdz—rdr; next, let the latitude of the body (or the angle between r and 

the plane of x, y) be X ; then tan X=-, and 

§ dz — zdg=r*d\, - 2 = sec 2 k. 

Making these substitutions, the expression for V becomes 

V=^+j , *(2mr 2 (A+9(r))-^)*+j^(^ 2 -c 2 sec 2 X)*. 

The integration in the second term cannot be effected till the form of the function 
<p(r) is given : that of the third term may be more conveniently performed after the 
differentiations with respect to c and k, as in the next article. 
29. The remaining integrals* of the problem are (art. 19.) 

dV dV a dV . , 

Performing the operations indicated, and observing that 

d\ i . _,/ AsinA \ 

, C sec 2 A.<& i . ,/ ctanA \ 

and JVF- c WA = ^ Sm \7F^> 

* It would perhaps be better to use the term " integral equations" here, in order to reserve the term " inte- 
gral" for the case of an equation involving only one arbitrary constant (see art. 23.). The equations — =a, &c. 

dk 

become "integrals" in this sense, when for k, c, and h, on the left, are substituted the functions of the varia- 
bles to which they are respectively equal (from (i.), (ii.), (iii.)). An *' integral" in this limited meaning is 
what is commonly called a " first integral," when the problem is considered as the solution of n differential equa- 
tions of the second order. And any equation obtained by combining " integrals " so as to eliminate a set of n 
of the variables x u a? a , ... x n , y t , y 2 , ... y n , of which no two are conjugate, corresponds to what is commonly 
called a " final integral." 



h 
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we obtain for the final integrals, 

m\rdr{2mr 2 (h-\-<p(r))— k i }-*=t-{-r (iv.) 

. . . / c tan A \ , . 

*- sm Atf^H (v - } 

-^-•{2mr^+9(r))-^}-*+sin-(^^)= a . . (vi.) 
Let j= cos / ; then == = cot/, and the equation (v.) becomes 

/c V ft ~~" c 

tanX= tan<.sin(^— 18), (v.a) 

which expresses that the orbit is in a plane whose inclination to the plane of*, y is /. 
Also (3 is evidently the longitude of the node, reckoned from the axis of x. 
The last term on the left of (vi.) becomes 

'sin A^ 



sin j 
Now if &be the "argument of latitude" or the angle between the node and the radius 

vector r, we have evidently sin ^=^^, so that the above term is simply 3-, and the 

am 

integral (vi.) becomes 

^_ a= ^j 2 mr 2 (A+p(r))-& 3 j - * (vi.a) 

30. To apply the above expressions to the case of the undisturbed motion of a 

planet, we have only to put <p(r)=^, where m is now the mass of the planet, and p 

the sum of the masses of the sun and planet, the origin of coordinates being placed 
at the sun. It would be useless to give the well-known expressions to which the in- 
tegrations now lead, my object being merely to obtain a set of normal elements. Now 
in this case we have (by well-known theorems), if a be the semiaxis major, e the 
excentricity, and / as before the inclination, 

h=^, k= ^a(l-e 3 ), 



and therefore c= Vpa(\— e 2 ).cos/. 

Also, if we take for the inferior limit of the integrations in (iv.) and (vi.a) the 
minimum value of r, or the perihelion distance, it is plain that a will be the longitude 
of the node, reckoned from the perihelion in the plane of the orbit, and — r the time 
of perihelion passage. Thus we have the following six elements, arranged in con- 
jugate pairs: — 

—-^5 —(time of perihelion passage) 

lid 



\/fj!,a(l — e*), (angle between node and perihelion) 



^pa{\— e 2 ).cos<, (longitude of node). 

o2 
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It is obvious that we may change the signs of the first pair. And generally, that 

iff, g be any two conjugate elements, we may substitute for them kf, -g, where \ is 

any determinate constant, i. e. not a function of the elements*. 

The above elements coincide with those given by Jacobi. My object has been 
merely to illustrate a mode of obtaining them which seems capable of useful applica- 
tions. 

31. As a second example I shall apply the method to the case of the motion of a 
solid body about a fixed point. 

Let the fixed point be taken for the origin, and the principal axes of the body 
through that point for the axes of*, y, z. Let |, q, £ refer to the same origin and to 
axes fixed in space ; a, b, c being the direction-cosines of the axis of x referred to 
the fixed axes of £, n, £, and d, b', d ; a", b", c" being respectively the direction-cosines 
of the axes of y and %. Let 6 be the inclination of the plane of x, y (or " equator") 
to that of |, n (or " ecliptic") ; 4- the longitude of the node, reckoned from the axis of 
|, and <p the right-ascension of the axis of x. Then if A, B, C be the Moments of 
Inertia, and p, q, r the angular velocities, about the axes of x, y, z, the expression for 

the vis viva is T=-x(Ap 3 -^-Bq 2 -\-Cr s ), where 

p — — & cos <p—>p' sin <p sin & 
q=tf sin Q—ty cos <p sin 

rrrrp'-f'^'cOS^. 

Let u, v, w be the variables conjugate respectively to 6, <p, 4>, so that 

dT dT dT 

dd' df' dty 

the following expressions will be found without difficulty : 

. , sin <p , „ 

Ap=—u cos p -}--=— f (v cos d—w) 

„ . , COS , 

iiq=u sin P+^-jj (» cos 6—w) 
Cr=«. 

Considering at present only the case in which no forces act, we have the integral of 
vis viva T=h, which becomes 

j\-u cos <P+-^j (v cos 6-w)j 

. 1 / . . COS a \ a 

+ 5 (usmp+-^(vcosd-w)) 

+^v 2 =2h. : (i.) 

* More generally, we may substitute for/, g any two functions of them,p, q, such that 

dp dq_dp dq _, 

dg df df dg~ 
a condition which requires the solution of a linear partial differential equation for the determination of one 
function, if the other be assumed. But on the subject of the transformation of elements see below, arts. 34, 35. 



DIFFERENTIAL EQUATIONS OF DYNAMICS, ETC. 101 

The three integrals which express the conservation of areas, namely, 

Aap •+- Ba'q + Ca"r = e, 

Abp+Bb'q+Cb"r=f, 

Acp+Bdq+Cd'r=g, 
become, after simple reductions, 

, sin ij/ . 

— u cos ip —•■ ^— r (v — w cos 6) = e 

. cos ty , „. ,, 

— »sin^/+- s ^y(t>— wcos^)=/ 

Let e 2 +/ 2 +5 ,2 =^ 2 ; we have, adding the squares of these three equations, 

2 1 2 1 ( V — WCOS0) 2 , 

" '+«>'+■ sin gfl =fi (11.) 

and we may take the three equations (i.), (ii.), and 

™=9 (iii.) 

as three normal integrals; the conditions 

[>,A]=0, [h,Jc]=0, [A,sr]=0 
being obviously satisfied. 

These three equations determine u, v, w as functions of 0, <p, ip ; and supposing the 
three former variables to be explicitly expressed in terms of the latter, we should 

have at once the three partial differential coefficients -^ — > -^ ; the determination 

r dd dtp aty 

of V would therefore depend upon simple integration, and the remaining integrals 

would be given by means of the three equations 

dj_ dV_ dj_ 

dh— t+T > dff— Ci > dk— c " 

r, c„ c 2 being new arbitrary constants. 

In the general case, however, the algebraical solution of the equations (i.), (ii.), (iii.) 
is impracticable, since the elimination of v and w leads to an equation of the fourth 
degree in u ; nor does it seem possible to evade the difficulty by choosing a different 
combination of integrals, since it may be shown that the necessary conditions cannot 
be satisfied unless two at least of the combinations chosen are of the second degree 
in u, v, w. 

32. Mr. Cayley has given* a solution of this problem, which, though differing 
totally in form and method from the above, resembles it in arriving exactly at a cor- 
responding point. For in Mr. Cayley's equations (27.), (28.), <X> and y are to be 
expressed as functions of v ; but this requires the algebraical solution of the system 
(18.) for p, q, r, and is therefore impracticable. (The two equations (i.), (ii.) of the 
* Cambridge and Dublin Mathematical Journal, vol. i. p. 167. 
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last article are merely transformations of the two first of Mr. Cayley's (18.); and (Hi.), 
though not identical with the third, is of the same degree ; so that the algebraical 
difficulty is precisely the same in both methods.) 

33. If we suppose A=B, the algebraical difficulty disappears, and the solution of 
the problem can be explicitly completed. But on account of the importance and 
interest of this case I shall make it the subject of a separate section, in which it will 
also be shown that the solution of the general case may be made to depend upon it, 
by means of the variation of elements. (See Section III.) 

34. Suppose any complete normal solution of the system of differential equations (I.), 
art. 14, be known, i. e. a solution involving the 2n elements 

a„a t , ... a n , b 1 ,b SS) ...b n 

which satisfy the conditions (23.), art. 9; then an infinite number of other sets of 
normal elements can always be found. 

For if we determine the 2» quantities «„ ... a„,ft, ... ft,, as functions of o„ &c., i„ &c. 
by the 2n equations 

dA t dA n 

dJ ( ~ * d^-P* 
where A is any arbitrary function of 

o v , a 2 , ... a„, « 15 a 2 , ... a B , 

it is obvious that the whole of the reasoning by which the formulae (19.), art. 7, were 
established may be repeated, merely putting A in place of X, and a, (3 instead of x,y. 
And repeating in like manner the reasoning of art. 9, mutatis mutandis, it will follow 
that if/, g represent any two of the 2» quantities a 1} &c, ft, &c, the expression 

v dtf, 9 ) 

will be equal to unity if/, g be a pair of the form «,, ft, and will vanish in every other 
case. But it was also shown ((25.) art. 9) that the above expression is equivalent to 
— [/, g] ; it follows then that 

[*j, ft] = - 1 , [«„ «,] = [«,, ft] = [ft, ft] =0* ; 

or, in other words, that a„ ... a„, ft, ... ft„ are a new set of normal elements. 

This method however can hardly be of much use in practice, because we cannot 
(at least without the solution of partial differential equations) determine what form 

;* I shall have occasion to refer afterwards to M. Desboves' Memoir in Liouville's Journal, vol. xiii., 
" Demonstration de deux theoremes de M. Jacobi." But it may be observed here that the proposition in the 
text is not the same as that expressed by the same notation in the memoir alluded to, p. 400. For M. Des- 
boves uses the symbols [sti, «,] in a different sense. His theorem, in the notation of the present paper, is 

^§^=1, or =0, 
according as/, g are of the form a,-, fy or not, which is easily established without the help of relations analogous 
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of the function A will cause any of the new elements to be given functions of the old. 
But the problems most likely to occur may be solved in another way, as follows. 

35. Assuming for the set a„ a 2 , ... a„, given functions of the set a„ a 2 , ... a n only, 
it is required to find |3„ ... (5 n . 

(It will be observed that the conditions [a u «,]=0 are necessarily satisfied in this 
case by virtue of (25.), art. 9, since «„ &c. do not involve b x , &c.) 

It is plain, that if the principal function X had been found from the n integrals 
a„ a 2 , ... a n (as in art. 14.), it would be changed into that which would be found from 
the n integrals a„ a 2 , ... a M merely by introducing the expressions for «„ ... a n in 
terms of a„ ... a n ; which expressions would be found by algebraical inversion of the 
assumed equations which give the latter set as functions of the former. Let X repre- 
sent the function X thus transformed ; we have then 

n dX dX da, dX da^ 

"' dtt { da 1 dai ' da 2 da-, 

«*.£+«£+■••+«£■ • • • • < 38 -> 

Thus (3 4 is determined as a function of the old elements, since ^, &c. may be ex- 
pressed in terms of the latter. In like manner we should have a set of inverse equa- 
tions 

which may be used instead of (38.). 

It is apparent that |3„ &c. will involve in general the elements a„ &c. as well as 

fe„ &c. 

Conversely, if we assumed for |3„ &c. given functions of the set b u ... b„ alone, we 

to (19.)> but would not answer our present purpose. I regret to use symbols with a meaning different from 
that which custom has to some extent sanctioned ; but there seemed to be only a choice of difficulties. 

Mr. SpottIswoode has suggested to me the employment of the symbols (analogous to Mr. Sylvester's 

" umbral" notation) 

u, v,u> 

d d d 
_ dx dy' dz* 

instead of those which I have used, namely, 



r^, y, z, ...~j 
\_d lt d 2 > d t ..._i 



d(x, y,z, . . .) 

If these were adopted, the two forms (p, q), [p, q] might be used without confusion in their usual significations. 
See note to art. 9. But although the " umbral" forms are more suggestive of the properties which belong to 
the above expressions as determinants, the other forms bring more into view the analogies which connect them 
with the differential calculus ; and therefore, for the purposes of this paper, I have preferred them. And it is 
perhaps better, for the present, that different notations should be tried, than that any attempt should be made 
to fix upon a definitive system for subjects so recent as those connected with the theory of determinants. 
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should have, for determining a„ &c, either of the systems 

a< =2,(«,!), *=*(<*g) (40.) 

We might obtain in this way an indefinite variety of sets of elements for the case 
of elliptic motion, beginning with those given at the end of art. 30. But it will be 
better to defer this illustration till after the discussion of the Method of the Variation 
of Elements, which will form the subject of a future Section. 

36. It results, from the investigations of this and the preceding Sections, that if a 

set of n integrals a„ a 2 , ... a, be given, satisfying the -^—z — - conditions [a iy aj=0, 

the determination of n more integrals b x , ... b n , constituting, with the given ones, a 
complete normal set, is a determinate problem, admitting of a unique solution, and 
always reducible (setting aside algebraical difficulties) to quadratures. 

But if, out of a complete normal set, n be given of which one or more pairs are 
conjugate, then the completion of the set is no longer a determinate problem, since 
the remaining n integrals, containing also one or more conjugate pairs, admit, to 
some extent, of arbitrary transpositions and combinations, as is evident from con- 
siderations similar to those employed in arts. 13 and 35. Hence we should expect 
a priori that the problem would require the solution of partial differential equations. 

It appears, indeed, at first sight, that having any n of the elements given functions 
of the variables, the relations established in art. 9, with the others included in the 
formula (21.), art. 9, would furnish more than a sufficient number of equations to 
determine explicitly all the partial differential coefficients of the remaining elements 
in terms of the variables*, at least in the case in which the principle of vis viva sub- 
sists, and the given integrals do not contain t. But it is certain from the above con- 
siderations that this cannot be the case, and therefore that the equations furnished by 
those conditions cannot be all independent. I have not at present attempted to show 
this directly, though it would probably be easy to do so. 

Note on art. 2, Section I. 

The theorem established in this article may be more shortly demonstrated as fol- 
lows : — 

Since dX(x^ t ) = %{x t dy t ) + 2 ( (y Ai) 

* The conditions [«,-, ij] = l, [«,-, fy]xs0, [*, bj]=0, [au a?]=0 will give, as is easily seen, n ( n ~ ' -j-n* 

equations ; and the analogous conditions (21.), art. 9, in which the summation refers to the numerators of the 
differential coefficients, will give the same number, so that upon the whole we shall apparently have 3b 2 — a 
equations, to determine the 2» ! partial coefficients required. 

It is not difficult to make mistakes in this subject. I was for some time under the impression that the pro- 
blem could be solved when any n independent integrals were given. Even the illustrious Jacobi himself 
appears to have been misled, at first sight, as to the consequences of Poisson's theorem (art. 22.). See the 
beginning of M. Bertband's Memoir mentioned above ; I do not know the fact from any other source. 
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and 2 i (y i (ki)=dX (by (5.)), 

we have 2,(avfy,)=d(— X+2,(a?«y,)), 

an equation which must become identical when a?„ x 2 , &c. on each side are expressed 

in terms of y n y 2 , &c. But the right side being then a complete differential of a 

ft f> j fif ' 

function of y„ y 2 , &c, the left side must be so also ; hence the conditions j- =T i must 

subsist. The investigation of art. 2 shows that they do subsist, and is therefore 
perhaps to be preferred. 

SkCtion III.— On the Equations of Rotatory Motion. 

37. In this supplementary section I propose further to exemplify the preceding 
theory by exhibiting the application of it to the problem of rotation in a more detailed 
form than was consistent with the plan of the former part of this essay. For this 
purpose it will first be desirable to anticipate the subject of a future section, so far as 
to give a concise deduction of the method of the variation of elements in its simplest 
form. 

38. Variation of Elements. — Suppose a complete normal solution of the system of 
differential equations 

*=*■ * ; +S=° • • • ■ w 

has been obtained, so that we have 2n elements, divided into two conjugate sets 

ffllj #2, ••• O-n 5 "15 "25 •'• "n 

as in the former articles, so that 

[a,, ft*] = 1 , [a,, a,] = [b„ ft,] = [a,, ftj =0. 

It is required to express the solution of the system 

,_dZ,dn ,,dz<m„ 

in the same form by means of variable elements. The disturbing function Cl may be 
a function of all the variables x H &c, y u &c, and may also contain t explicitly. 
In the undisturbed problem we have «-=0, ft,=0 ; i. e. the equations 

5+LZ.«a«=0,f+[Z,*a=0 (e.) 

(see art. 22.) subsist identically when x 1} &c, y u &c. are expressed in terms of the 
elements and t. 

In the disturbed problem, x 1} &c, y 1} &c. are to be the same functions of the 
elements and t as before ; hence the equations (e.) continue to subsist identically, and 
therefore the values of d i} ft-, namely, 

ft;=§+[z/&j+[n,*a, 

MDCCCL1V. P 
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become simply a' { —[Cl, a,], &•— [O, &,]. 

In these expressions O, a 4 , 6 f are supposed to be expressed in terms of the variables. 
Now 

but, by equation (26.), art. 10, this is equivalent to 

d{D,, a t ) 
- Z Sd(bj,aj)> 

in which fi is expressed as a function of the elements and t ; and this last expression 
obviously reduces itself to the single term — — • In like manner the expression for 

[CI, bi] reduces itself to +-r— ; thus the equations for determining the variation of 
the elements are 

a ^~^ *<=?? • • • • < R ) 

in which CI is to be expressed as a function of the elements and t*. This will be a 
sufficient account of the method for our immediate purpose. 

39. The following propositions in spherical trigonometry will be required. If a, b, 
c be the sides, and a, /3, y the opposite angles of any spherical triangle, then 

(cos a. + cos B)' z , 

; — — 1 — cos a cos p 

co S (a+b)= l - cos r (40.) 

smasin/3 v ' 

_ ( cos«-^os^ +1 _ cosaco 

COS (a-b)= 1+C08 r . g ; (41.) 

' sin a sin jS v ' 

and if the sides be considered as functions of the angles, then 

da „„„ db i „odc /.„\ 

_=co S y-+cos^^ (42.) 

da „ db , r> dc . . „ > 

-^cosy-^+cos/3^- (43.) 



T 



The two last are easily verified ; but as the others are not so obvious, I shall give the 
demonstration. Putting x for the expression on the right of the equation (40.), we 

* The history of these remarkable formulae may, I believe, be stated as follows. They were first discovered 
by Lagrange in the case in which ai, bs were the initial values of #,-, yi, and £1 contained w u &c. but not y,, &c. 
They were extended by Sir W. R. Hamilton to the case in which Q, contains both sets of variables ; and 
finally, by Jacobi, to the case in which a v &c, b lt &c. are any system of conjugate elements. Jacobi how- 
ever does not appear to have published a demonstration of them, and the only one which I have seen is by 
M. Desboves, Liouville's Journal, vol. xiii. p. 397, and differs essentially from that given in the text. 
Sir W. R. Hamilton has pointed out the circumstance, that when Q, contains both sets of variables, the vary- 
ing elements determined by the formula (E.) are not osculating. 
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easily obtain 



co8»5=£ sin^-cos*^ 



2 



cos 2 — — - cos 2 — — - — sin 2 ! - 



-*—$ a + fi + y a+/3— y 
cos 2 — - cos - — - COS - - 



2 « + /3 a— (3 + Y (3+y — a 
COS 2 —— 1 - COS ~ ^COS- fr 

Z Z 1 



COS 2 — — - 



- tan 2 -x= tan 



2 a + i 



"cos 2 ^ 2 2 



whence it is plain that a?= cos (a+6), and in like manner may the equation (41.) be 
established. 

40. Returning now to the problem of rotation, and supposing, for convenience, 
that the question refers to the motion of the earth about its centre of gravity, the 
following will be the signification of the symbols employed. 

A, B, C are the moments of inertia about the principal axes of the earth, viz. the 
axes of x, y, z ; the last being the polar axis, and the arrangement being such that 
the positive direction of z is to the north pole, and that the positive axis of x follows 
that of y in the actual rotation about the polar axis : p, q,r being the angular velo- 
cities about the three principal axes, the usual convention will be adopted as to their 
signs ; so that in the actual case r is positive. The arrangement of the fixed axes of 
I, n, Z, is supposed similar to that of x, y, z, the plane of g, n being a fixed ecliptic, and 
the axis of g the origin of longitudes unless another origin be expressly indicated. 

Then 6 is the oliquity, 4» the longitude of the vernal equinox, and <p the right 
ascension of the axis of x ; all referring to the fixed ecliptic. 

Let the "principal plane" signify that which, in the undisturbed problem, is the 
"invariable plane." Then i is the inclination of the principal plane to the fixed 
ecliptic, and j is the inclination of the equator to the principal plane. 

In the case of the earth, A is nearly equal to B, 6 never differs sensibly from i, and 

j is therefore always small. But these conditions are not supposed in what follows. 

It is assumed however that C is the greatest of the three moments of inertia. These 

conventions, in which it is very desirable to avoid any ambiguity, may be illustrated 

by the annexed figure, in which O represents the origin of longitudes. 

The angles of the spherical triangle formed by the intersection of the three planes 
with a spherical surface are i,j, v—d\ and the sides opposite to them will be denoted 
by I, J, 6. Thus we shall have 

cos I— C0Si '- c0S ./ C0S fl , cos j^ cosy-coszcosfl 
siry'sinA sin i sin 8 

cos 6=. cos i cosj— sin i sin; cos 0. 
p2 
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And, in the figure, OT=i£, and p is measured from T in the direction indicated by 




the arrow, which is also the direction of the rotation about the polar axis. Moreover, 
if the direction-cosines of the axes of x, y, z referred to the fixed axes, be respectively 
a, b, c; a', V, d ; a", h", c", we shall have 

a = cos -^ cos <p<— sin 4> sin p cos 6 

a! = — cos 4 1 sin p— sin 4* cos p cos 

a"=z — sin 4- sin 

b — sin ^ cos <p+ cos ^ sin p cos 6 

b' = — sin ip sin <p+ cos ^ cos p cos 8 

b"=s: cos 4 sin # 

c =— smp sin 4 

c' = — cos p sin 

c"=cos0 

P=—0' cos p— •>£' sin p sin & 

§-=# sin <p— 4' cos p sin 4 

r=p'+-4/ cos 6, 

hence we obtain the expressions for u, v, w employed in art. 31, viz. 
«=— ,= — Ap cos p-\-Bq sin p 

dh /-, 



w=^j- = —Ap sin p sin 6— Bq cos p sin 4+0 cos 4, 
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from which the following also are easily deduced : 

u=s-j-6'-\—^—(0 t cos2p+^'s\n&sm2p) 
w=?C(p'+%|/cos<9 

A I JJ ^ g 

w=i— g— 4/ sin 2 0+C cos (<p'+4' cos 0)H — g— sin 6 (0 sin 2<p — %^' sin & cos 2<p). 

41. Resuming the three integrals (i.), (ii.), (iii.) of art. 31, we may put the first in 
the following form : — 

Z+0=A, ,. (i.) 

in which 2Z=2( A+l) (*+ sin2fl ; )+c 

and the other two are, as before, 

M> 2 +M a + sin ag =^ 3 (ll.) 

">=#, (iii.) 

in which & is the sum of areas on the invariable plane, and g the sum on fixed ecliptic ; 

moreover vz=Cr is the sum of areas projected on the plane of the equator; hence 

we have 

gz=kcosi, vz=k cosj\ 

It has been seen that the complete solution of the problem is impracticable in the 
general case, on account of an algebraical difficulty. If however we suppose B=A, 
this difficulty disappears ; and after completing the solution on this supposition we 
may take account of the terms arising from the inequality of A and B, by treating 
the function denoted above by O (equation (i.)) as a disturbing function, and apply- 
ing the method explained in art. 38. Thus when the action of disturbing forces is 
considered, the whole disturbing function will consist of two parts ; one depending 
upon the forces, and the other the function which has just been assigned, and of 
which the effect, as will be seen, is extremely simple. 

42. We proceed then first to complete the solution on the supposition A=B. The 
three integrals (i.), (ii.), (iii.) give in this case 

1> 2 ~C3£(& 2 — 2AA), w—g ........... . (44.) 



M= sn7&j k 2 —v*~w 2 +2vw cos 6—W cos 2 ffV 



in which latter expression the above constant values of v and w are to be introduced. 
We may put, as before, g=# cos i, v=sk cosj,j being now constant ; and the expres- 
sion for u becomes 

w== sm0I 1 — cos2 *"" cos2 i+ 2 c °s * cosj cos 0— cos 2 #l . . . (45.) 
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and we shall have (art. 19.) 

V=k(yp cos i-\-<p cos j)-\-yud6; . (46.) 

and we will take 

h, cos i, cosj 

for normal elements*, so that k is to be considered as a function of these elements, 
given by the equation (see (44.)) 

hi _ 2ACfl U7 v 

ft -C-(C-A)cosV K '' } 

It is to be observed, that, according 1 to the hypotheses admitted above, k is positive ; 
also, the expression for u at the end of art. 40 becomes, in the case now considered, 
u=Ad'. Thus u has the same sign as 6 ; and since 6 is evidently comprised between 
i—j and i-\-j, if we suppose i andj both acute (as in the figure), sin 6 is always posi- 
tive ; hence in the expression (45.) for u, we have to attribute the sign -j- or — to 
the radical, according as 6 is increasing or diminishing, or according as © is between 
o and t, or not ; thus, in the position represented in the figure, the negative sign 
must be taken. 

43. If we put +Q for the radical in question, the expression for udd is easily trans- 
formed into the following, namely, 

,. . ksinidO f 1 (cosj— cosi) 2 1 (cos^'H- cosi)' 

ud6—± q— |1— g 1-cosS 2 l+cos0 

in which it is evident that the part within brackets is positive upon the whole, but 
each of the two last terms is essentially negative. The integration is now easily per- 
formed, and the result is 



J* 



\ud6=±kV, 

where the sign is that which belongs to the radical Q, and P is given by the equation 

~, . cos 6 — cos i cos j 

P= COS" 1 - =—=-= ; ~ 

sin i siry 

(cos/— cos/) 2 , 

.i — ^l _; — 1+ cos « cos; 

1 , ., , 1 — cos J 

— ~(C0S7 — COS*) COS" 1 — . . . ■ 

2 V J > sin « sin; 

(cos ;'+ cosi) 2 , 

- — ~ j- 1 — 1— cose cos; 

. 1 , . . ., _. 1 + cos J 

+h(cos?+ cos^) cos ' • • • ■■ ■ , 

1 2 V J ' ' sin 2 sin; 

+ an arbitrary function of i,j, h. 
This apparently complicated expression has a very simple geometrical signification ; 

* Since*?, h, k are normal elements (i.e. satisfy the conditions [g, h] =0, [h, k]=0, [k, g]=0, art. 31.), any 
three independent combinations of them are obviously normal also. 
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for, referring to the figure, and using the theorems (40.), we see that it is equivalent to 
P= cos-H- cos Q)+ C0S ^ g C0S * cos-'(cos(I+J)) 



COS J— COS J 



cos- 1 (-cos(I-J))+K, 



where K is put for the arbitrary function. Now the expression for udd (from which 
this is derived) shows that the three terms in the above value of P must be so inter- 
preted that the differential coefficient of the first (with respect to 8) shall be positive, 
and those of the two others negative. These conditions will be satisfied by taking* 

+P = ,_0 + £2!4^i(I +J) 

- coaj - co %-{l-J))+K 

(in which the upper sign is to be taken when is between o and or, and the under 
sign when is >w). 

Hence, assuming the arbitrary K so as to destroy the constant part of the expres- 
sion, we have, without ambiguity, for all values of the variables, 



s- 



\udd=k(®— I cosj— J cos i), 

so that, finally, 

V=/c{(4/-I)cosj + (?)-J)cosi+0} (48.) 

It will be observed that without attention to the proper interpretation of ambiguous 
symbols, a completely erroneous expression for V might have been obtained. 
44. The final equations will be (art. 19.) 

dV., JI dY _ 

dh — t-\-r, dcosi —<*, dcosj—P' 

t, a, (3 being three new arbitrary constants, namely, the elements conjugate respect- 
ively to h, cos i, cosj. 

In performing the differentiations, it is to be remembered that I, J, do not con- 
tain h; and that, by the equations (42.), (43.), art. 39, the terms arising from the 
differentiation of I, J, with respect to i and j, disappear identically, so that these 
functions may be considered as exempt from differentiation. Also we have 

dk k dk 

WiT" 2A' dcosi~ ' 

dk _ (C-A)k s cosj 
dcosj~~ 2ACA 

* In the figure, as 8 diminishes (i and,/ remaining constant) 9 increases, 1+ J increases, and I— J increases or 

• i—j 
t_j e sul— o 
diminishes according as,/ i, since tan — — - = tan— __. 

sin^-^- 
2 
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(see equation (47.)), and the final equations become, after simple reductions, 

«cos«'+/3cos,/ . k 

■ r + 5 (*+r) 



e=s- 



^-I=f-^-gJAC0Sj.(i+r) 



^-J=T 



(R.) 



These equations comprise a normal solution of the problem. The first gives imme- 
diately 



(k & 



COsi-t-ffCOsA 

k ) 



(see art. 40.) ; and since I, J are given explicit functions of 6, the three variables 
6, <p, 4> are determined explicitly as functions of t. The third equation (R.) simply 
expresses that the invariable plane intersects the ecliptic in a fixed line, whose longi- 



tude is t- 



45. Let us now introduce the supposition that A and B are unequal, and that the 
body is acted on by disturbing forces. 

We must (see art. 41.) putovA-Hj) mstea ^ °f T * n tne e q uatlon s (R.) of the last 

article ; these equations will express the solution of the problem, the elements being 
now variable, and determined as functions of t by the system of equations 






d® 



l>-r 



d® 



dcosi r dcosj 

where 3> is the disturbing function, expressed in terms of the elements and t* 

46. If there are no disturbing forces, O reduces itself simply toO (art. 41.), which 

is now to be transformed by means of the equations (R.), art. 44, as follows. 
Since v=&cosj, and w=kcosi, we have 

pcosfl — w , cos i— cos; cos 5 , . . T 

sin i =-* STS =-.* S1 ^ C0S L 

Also the expression for u, art. 42, is easily put in the following form s 
u= -r— J sin 2 i sin 2 j— (cos 6— cos i cosj) 2 

k 
~ ""sinfl Sm * Sm «^' Sm ® 

(with respect to the sign, see art. 42.). And since 

sin® sin I 

sin 6 sin?' 
uss— k sin j sin I, 



this becomes 
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Introducing these expressions in the value of ft (art. 41.), we find 

n= -f(i-B) sin2 i cos2 ^- I ) ; ( 49 -) 

and when <p— I is expressed in terms of the elements and t (see equations (R.), art. 44, 
in which ^l a+b) * s now to ^ e Wl 'i tten f° r t) 5 this becomes, finally, 

n =-i(I-s) shl2 i-eos2[f-(i(i+i)-i)/ccosi.(^+r)]. . . (ft.) 

47- The above expression for ft does not contain the elements i, a, ; hence, when 
there are no disturbing forces, we shall have (cos?')'=0, a'=0, or i and a are con- 
stant ; also 



,, dk d£l dk dD, 

dh dr d cosj dfi ' 



dk dk 



an expression which is easily found to vanish identically (see the values of , 

in art. 44, observing to put -A ^+g ) fc> r r ) • Thus k is also constant ; and the " prin- 
cipal plane" is still the " invariable plane," as we know a priori. 

48. If we now suppose the attraction of another body to be introduced as a disturb- 
ing force, we shall have to take for the disturbing function 

<D=ft-P, 
where ft is the same as above, and P is the potential of one body upon the other, ex- 
pressed as a function of the elements and the time*. And it follows from the remarks 
of the last article, that the variation in the position of the principal plane depends 
wholly upon P, and not upon ft. 

I shall here conclude this part of the subject, as it would be beyond the scope of 
this essay to enter into the details of any of the various problems which might be 
taken in illustration of the theory, such as those which relate to precession and nuta- 
tion, or to the motion of the moon about its centre of gravity. The investigations of 
this section have been introduced, because the results, so far as they go, appeared 
interesting in themselves, and afforded a remarkable example of the application of 
the general method. 

P.S. Since the last sheets of this essay were in type, I have seen for the first time 
two papers by Professor Brioschi, in Tortolini's Annali for August and October 
1853, of which the titles are 'i Sulla variazione delle costanti arbitrarie nei problemi 
della Dinamica," and " Intorno ad un teorema de Meccanica." I have not had an 
opportunity of examining them sufficiently to judge how far any of the preceding 
investigations may have been anticipated in them. 

June 7. 

* The variables which determine the position of the disturbing body are supposed to be given explicit func- 
tions of t. 
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